Pade-Improvement of CP-odd Higgs Decay Rate into Two Gluons 
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' We present an asymptotic Pade-approximant estimate for the four-loop coefficients 

within the hnear combination of correlators entering the recently calculated decay rate 
of a CP-odd Higgs boson, with an assumed mass niA = 100 GeV, into two gluons. 
All but one of these coefficients are shown to be determined for arbitrary from 
the known 3-loop-order rate by renormalization group methods. Asymptotic Pade- 
approximant estimates for these coefficients are all seen to be within 12% of their 
correct values. The four-loop term in the decay rate for tjia = 100 GeV is estimated 
to be only 4.3% of its leading one-loop contribution. 
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The decay rate into two gluons [g) of a CP-odd Higgs boson [A) occurring within a 
two-Higgs-doublet extension of the standard model has been calculated to three loop order 
by Chetyrkin, Kniehl, Steinhauser, and Bardeen Their result is expressed in terms of a 
linear combination of the imaginary parts of three different correlators: 
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For Nc = 3 and n light flavours, the terms within the linear combination (2) are shown IQ 
to be [x = as/iT, L = ln(/i^/g^), L' = ln(m^/g^)] 
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One can combine these results to obtain the following series for the linear combination of 
correlators defined by (2): 
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The terms listed above arise entirely from the first two terms of (2), as the final term 
(C'|/m([02]^)) is 0{x^). The four-loop [(9(x^)] coefficients Cq — C3 in (8) are as yet unde- 
termined. All but Co of these can be obtained via renormalization-group (RG) methods. 
RG-invariance of the physical decay rate (1) and, consequently, the linear combination of 
correlators (2) implies that the function S[x,L,L'] in (8a) satisfies 
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S[x, L, L'] = 0, 



where 



p(x) = -(3oX^ - pix^ - p2X^ . . . 
Imt {x) = ~loX - 7lX^ - 720;^ . . . 
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If rrit is a pole-mass independent of the renormalization scale /z, then = 0. However, if 
rrit is a /i-dependent running quark mass, then 70 = 1, and subsequent 7j's in (11) are as 
determined in ref. P|. 

Substitution of (8a), (10), and (11) into (9) yields the following set of equations for the 
aggregate coefficient of x^^L^'^ to vanish: 



x^: ai - 2po = 0, 
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x^: hi - 3(3oao - 2/5i = 0, (13) 

x^L: 262 - 3ai/?o = 0, (14) 

x^ ci - AboPo - 3ao/3i - 2p2 + 2n-fo = 0, (15) 

x^L: 2c2 - 46i/3o - 3ai/3i = 0, (16) 

x^L^: 3c3 - 462/5o = 0. (17) 
The coefficients P0-2 in (10) for n light flavours are given by 

A = ^ - ^. (18) 
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the coefficients 71, 72, • • • in (11) do not enter (9) until 0{x^\ Using eqs. (8.b,c,e,f), (18), and 
(19), we see that eqs. (12-14) are explicitly upheld, thereby conffiming the RG invariance of 
(8a). The unknown coefficients Ci, C2, and C3 are obtained via equations (15-17): 

ci = 4822.88 - n(llL' + 884.455 + 270) + i^U + 45.1091^ - 0.591921n^ (21) 

V 2 3/ V 8 48 12 y ' ^ ' 

1 /II , , 

For the physical case of n = 5, = 175.6 GeV [the t-quark pole mass (70 = 0)], with rriA 
chosen as in [1] to have a reference value of 100 GeV ^ we find that 

ci = 1411, C2 = 438.4, C3 = 28.16 . (24) 

The coefficient cq is RG- inaccessible to order x^ . 
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The four-loop correlation-function coefficients cq-s can be estimated using asymptotic 
Pade-approximant methods as delineated in ref. Q]. Given a correlation function of the 
form 

n(x) = F{x) [l + Rix + R2X^ + R^x^ + ...], (25) 

with only coefficients Ri and R2 known, the simplified asymptotic error formula (utilized in 
to estimate [3^ from /?o-2) 

d[^|1] d a 

characterizing the [A^|l] Pade-approximant prediction for Rn+2, yields the following predic- 
tion for 0: 

R^ = 2Rl/{Rl + R^R2). (27) 

Comparing eq. (25) to (8a), we see that the coefficients -Ri, i?2, -R3 are necessarily functions 
of L = ln(/iVg2). 

i?i = ao + aiL, (28a) 
i?2 = &o + biL + b2L', (286) 
^3 = Co + CiL + C2L^ + cgL^ (28c) 
Consequently, we can obtain co_3 from the moment integrals 

Nk = {k + 2) dww^+^R^{w), (29) 

JO 

where w = = -ln(w)].[| Exphcit substitution of (28c) into (29) yields [Q 

A^_i = Co + ci + 2c2 + 6C3, (30) 
113 

No = Co + -ci + -C2 + -C3, (31) 



^Such integrals characterize the 0{x^) contributions to the /c*'* finite-energy sum rule integral 
over the correlator (25), where in (29) corresponds to t, and where /i^ in (29) corresponds 
to the continuum threshold sq. 
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Numerical values of A^_i, A^o, -^i, and N2 can be obtained through explicit use of the Pade- 
motivated estimate (27) within the integrand of (29) with and R2 given by (28a) and 
(28b). Within these latter two equations, the coefficients Oq, Oi, ^i; ^2 are as given by 
(8.b-f). We choose n = 5 and L' = 2 ln(1.756) to facilitate comparison with the true RG 
values (24) for = 100 GeV , and find that 



The relative errors of the above Pade estimates for ci, C2, and C3 from their true values, as 
given in (24), are respectively -9.7%, +3.2%, and -11.3%. 

An alternative method for extracting co_3 is to fit Rsiw), as obtained from (27), to the 
form of (28c) via least-squares minimization of the following function: 



A^_i = 3310.4, A^o = 1863.8, Ni = 1512.6, N2 = 1359.2 . 



(34) 



We substitute these values into (30-33) to find that 



Co = 981.7, ci = 1274, C2 = 452.3, C3 = 24.96 . 
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The values of Co_3 which minimize are 



Co = 978.9, ci = 1285, C2 = 445.0, C3 = 26.03, 



(37) 
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in very close agreement with the values (35) extracted from the moment integrals (29). 
Moreover, is equal to only 4.7 at this minimum. The near cancellation of the 0(10'') lead 
term in (36) is indicative of the precision of the fit obtained between (27) and (28c), as is 
evident from Figure |l|. Relative errors of {ci, C2, C3} obtained from (36) with respect to their 
true RG-determined values (24) are respectively -8.9%, +1.5%, and -7.6%, confirming the 
usefulness of the asymptotic Fade approach in estimating four-loop order contributions to 
the correlation function (25). 

The iriA = 100 GeV estimate for cq can be improved somewhat by using the correct 
values (24) of ci_3 within equation (30), the lowest moment integral estimated in (34) by 
asymptotic Fade-approximant methods. We then find that 



Identically the same result is obtained by minimizing with respect to cq after explicit 
incorporation of the correct (RG) values (24) of Ci_3 into (36). There is a 15% discrepancy 
between (38) and the estimates for cq in (35) and (37), indicative of the magnitude of 
anticipated relative error with respect to the true value for cq. It is hoped that these estimates 
can be tested against an exact 4-loop calculation in the not-too-distant future. 

The 4-loop correction to the CF-odd Higgs decay into two gluons, as determined to 3-loop 
order in [0], is found from (1), (8) and (38): 



Given n = 5, = 100 GeV, a^^^(m^) = 0.116 Jll, and mt = 175.6 GeV, the square 
bracketed expression in (39) for successive-loop corrections is [1 + 0.680 + 0.226 + 0.043] . The 
first three numbers are as calculated in ref. IH; the final (underlined) term is obtained from 
the asymptotic Fade-approximant estimate for cq in (38). This estimate is further indicative 
of a progressive decrease in the ratio of successive terms in the A 2g decay rate, suggesting 
that if such a CF-odd Higgs were discovered, a perturbative calculation of its 2-gluon decay 



Co = 3310.4 - 1411 - 2(438.4) - 6(28.16) = 854. 
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rate could lead to a phenomenologically testable value [e.g. 1 + 0.680 + 0.226 + 0.043 + ... 
K, 2.0 for uiA = 100 GeV]. Such a Higgs characterizes the two-doublet version of electroweak 
symmetry breaking anticipated from supersymmetric extensions of the standard model, as 
first noted over two decades ago 0. 
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Figure 1: The ratio of the Fade prediction R^'^^{w) (27), and the x^-minimizing fitted form 
rI^*{w) (28c) as a function of w. 
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